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Abstract 

Boot  locus  analysis  techniques  are  used  to  design  an  active 
thrust  vector  control  system  to  provide  the  directional  stability 
for  an  F-lll  without  the  vertical  stabilizer.  A  linear  analysis 
of  the  lateral-directional  modes  is  performed  for  the  airc  aft  both 
■with  and  without  the  vertical  stabilizer.  The  aircraft  with  the 
vertical  stabilizer  is  used  as  a  baseline.  Computer  time  histories 
for  discrete  atmospheric  turbulence  and  a  covariance  analysis  for 
random  turbulence  are  used  for  the  evaluation.  It  is  found  that 
the  thrust  vector  control  produces  a  response  as  good  as  or  better 
than  the  baseline  aircraft. 

Requirements  for  the  thrust  vector  deflection  and  rate  of  de¬ 
flection  are  generated.  The  lowest  possible  rate  of  deflection  for 
acceptable  flying  qualities  is  shown  to  be  1  rad/sec. 


THRUST  VECTORING  TO  ELIMINATE  THE 


VERTICAL  STABILIZER 

I.  Introduction 

There  is  increasing  concern  ove-  the  vulnerability  of  modern 
jet  fighter  aircraft  to  detection  by  radar  and  optical  means.  In 
looking  at  ways  to  reduce  detection  there  has  been  much  work  on 
electronically  masking  an  aircraft  to  radar  ,  although  these  techniques 
do  little  to  reduce  the  optical  profile.  When  one  views  a  jet  fighter 
from  a  top  view  (or  bottom  view),  the  wings  constitute  the  largest 
part  of  the  aircraft's  cross  section.  Since  they  have  been  designed 
specifically  to  give  the  aircraft  its  required  lift  and  flight 
characteristics,  there  is  no  apparent  potential  for  improvement  of 
the  cross  section.  From  a  side  view  of  the  aircraft  it  is  seen  that, 
besides  the  fuselage,  the  cockpit  and  the  vertical  stabilizer  are 
areas  in  which  to  improve  the  cross  section.  Reduction  of  the  height 
of  the  cockpit  will  not  help  very  much  and  would  reduce  the  pilot's 
visibility.  Reduction  of  the  size  or  elimination  of  the  vertical 
stabilizer  and  rudder  would  greatly  improve  the  radar  cross  section 
and  optical  profile.  The  problem  is  that  the  vertical  stabilizer 
provides  the  aircraft's  directional  stability  and  therefore  its 
reduction  or  elimination  could  cause  the  aircraft  to  be  directionally 
unstable.  Therefore,  a  means  to  provide  the  directional  stabilility 
and  yaw  characteristics  that  the  vertical  stabilizer  and  rudder 
provide  is  needed.  Thrust  vectoring  of  the  jet  engine  nozzle  using 
an  active  control  system  could  be  the  means. 


Another  problem  with.  the  vertical  stabilizer  is  that  at  the  high 
angles  of  attack  modern  jet  fighter  aircraft  are  capable  of,  the 
stabilizer  may  be  blanked  to  the  free  stream  by  the  wing-body  combin¬ 
ation.  An  aircraft  with  an  active  thrust  vector  control  (TVC)  would 
not  be  as  susceptible  to  reduced  lateral-directional  stability  at 
high  angles  of  attack. 

Furthermore,  if  the  vertical  stabilizer  could  be  reduced  or 
eliminated  ,  both  a  reduction  in  drag  and  weight  would  be  realized. 

Of  course,  these  benefits  must  be  compared  with  the  increased  com¬ 
plexity  and  weight  and  the  reduction  in  longitudinal  thrust  associated 
with  thrust  vectoring. 

Purpose 

The  purpose  of  this  thesis  is  to  determine  the  static  and  dynamic 
requirements  of  lateral-directional  thrust  vectoring  (magnitude  and 
rate)  adeqxiate  to  replace  the  vertical  stabilizer  on  a  specific  air¬ 
craft.  In  other  words,  this  thesis  will  determine  what  is  needed  in 
an  active  TVC  to  provide  the  same  directional  stability  and  control 
now  given  by  the  vertical  stabilizer  and  rudder.  Also,  the  stability 
and  control  advantages  of  such  a  concept  are  to  be  indicated. 

Aircraft 

The  aircraft  utilized  in  this  study  is  the  F-lll  Transonic 
Aircraft  Technology  (TACT)  program  aircraft.  The  F-lll  is  a  super¬ 
sonic,  variable  wing  sweep,  twin  engine  fighter  aircraft.  A  physical 
description  of  the  F-lll  TACT  aircraft  is  in  Appendix  1.  The  TACT 
program  (Ref  1:3)  was  a  supercritical  airfoil  research  program  to 
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demonstrate  the  performance  advantages  of  a  supercritical  wing.  Since 
wind  tunnel  data  for  this  aircraft  in  both  a  tail-on*  and  tail-off 
configuration  was  available  (from  the  Flying  Qualities  Group  of  the 
Air  Force  Flight  Dynamic  Laboratory) , it  was  a  convenient  choice  for 
this  present  study. 

Scope 

Linear  analysis  and  design  methods  were  used  to  approach  the 
problem.  Since  it  was  initially  assumed  that  the  longitudinal  pertur¬ 
bations  would  be  small,  only  the  lateral-directional  aircraft  modes 
were  examined.  The  limitation  of  this  assumption  is  discussed  in 
Chapters  IV  and  V.  A  full  design  study  of  the  aircraft  covering 
its  complete  flight  envelope  was  not  practical  for  this  thesis. 

Instead  the  study  was  restricted  to  two-steady  level  cruise  conditions. 

The  analyses  and  design  were  performed  at  cruise  conditions  of 
Mach  0.8  at  1*0000  ft.  altitude  and  at  Mach  2.0  at  30000  ft.  Since 
the  first  is  a  low  thrust,  low  dynamic  pressure  condition,  it  is 
believed  to  be  the  a  worst  case  cruise  condition  with  respect  to  damping 
of  the  dutch  roll  mode  with  the  TVC.  On  the  other  hand,  the  second 
flight  condition  is  a  high  thrust,  high  dynamic  pressure  condition 
and  therefore  represents  an  easier  design  case.  Table  1  summarises 
these  two  flight  conditions.  The  stability  derivatives  for  the  two 
flight  cases  in  the  tail-on  and  tail-off  configurations  are  shown  in 
Appendix  2. 

The  system  was  evaluated  using  both  discrete  and  random  atmo¬ 
spheric  turbulence  simulation.  The  tail-on  aircraiu,  including  Its 

*Tail-on  refers  to  the  aircraft  with  the  vertical  stabilizer  and 
tail-off  refers  to  the  TVC  aircraft  without  the  vertical  sta¬ 
bilizer. 
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Table  1 


Flight  Conditions 


Altitude  ('f't) 


1*0000  30000 


Airspeed  (ft/sec) 

Dynamic  Pressure  (lb /ft2 ) 
Density  (sluGs/ft3) 


771* 

176 

-1» 

5.873^10 


1990 
176U 
8.907x10 


Flight  Path  angle  (deg) 


0 


Angle  at  attack  (deg) 
Weight  (lb) 

Wing  Sweep  (deg) 

Center  of  Gravity  ($  MAC) 
Required  Thrust  (lb) 

Ixx  (slug- ft2) 

Ixz  (slug-ft2) 

f  -*  4U  2  \ 

jlz z  \yj.u / 


9 

2 

68000 

68000 

26 

58 

32.5 

32.5 

6350 

1*1*980 

89500 

59000 

3500 

1*720 

1*56000 

1*1*2000 

yaw  stability  augmentation  system  (SAS),  was  used  as  a  basis  for  com¬ 
parison.  A  sensitivity  study  was  performed  to  determine  the  effects 
of  thrust  variation,  feedback  gain  and  TVC  actuator  time  constant 
variation.  The  TVC  actuator  requirements  to  provide  acceptable 
flying  qualities  were  specified.  No  attempt  was  made  to  model  any 
specific  two-dimensional  TVC  nozzle.  Only  the  deflection  of  the 
thrust  vector  is  specified.  In  its  undeflected  position  the  thrust 
vector  is  assumed  to  be  along  the  aircraft  body  roll  axis. 

Organization 

The  second  chapter  of  this  thesis  presents  the  lateral-directional 
small  perturbation  equations  of  motion  utilized  in  this  study. 

Chapter  III  describes  the  design  of  the  thrust  vector  controller. 
Listed  first  are  the  flying  qualities  requirements  which  the  design 
is  to  meet.  The  block  diagram  of  the  yaw  control  system  is  shown  and 
then  the  root  locus  analysis  of  the  controller  is  presented.  The  feed¬ 
back  gain  selections  are  provided  along  with  the  transfer  functions  of 
the  closed  loop  system  with  the  yaw  SAS.  Then  a  comparison  is  made 
with  the  tail-on  yaw  SAS  and  verification  that  the  flying  qualities 
requirements  are  satisfied  from  the  eigenvalues  of  the  characteristic 
equation. 

The  results  of  the  evaluation  of  the  yaw  SAS  to  three  types  of 
input  are  presented  in  Chapter  IV.  A  computer  program  was  developed  to 
integrate  the  aircraft  equations  of  motion  to  provide  time  histories 
for  a  pulse  rudder  pedal  input.  A  comparison  is  made  with  the  tail-on 
aircraft.  Next,  the  response  to  discrete  atmospheric  turbulence  is 
shown.  The  model  for  the  turbulence  is  a  1-COS  discrete  side  gust. 
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The  aircraft  response  to  this  gust  is  found  using  the  some  integration 
program  as  before.  Again,  a  comparison  -with  the  tail-on  configuration 
is  made.  The  effect  of  random  atmospheric  turbulence  is  next  examined 
using  the  Dryden  form  of  the  power  spectral  density  and  a  computer 
program  which  does  a  covariance  analysis.  Results  of  a  study  of  the 
effects  of  the  •‘'"'nation  of  thrust,  gains,  and  actuator  time  constant 
(each  independently)  on  the  response  to  atmospheric  turbulence  are 
presented . 

Chapter  V  then  discusses  the  effects  on  the  dutch  roll  mode  of 
applying  limits  to  the  thrust  vector  excursion  and  slew  rate  using 
the  discrete  atmospheric  turbulence  input.  The  results  show  the 
minimum  thrust  vector  excursion  and  rate  requirements  for  the 
discrete  gust  environment  to  satisfy  the  dutch  roll  requirements. 

Conclusions  and  recommendations  are  presented  in  the  final 
chapter  of  this  thesis . 


II.  Aircraft.  Eauations  of  Motion 


This  chapter  presents  the  lateral-directional  small  perturbation 
equations  of  motion  (Ref  2:6.U7)  written  in  the  body  axis  system. 
These  equations  were  used  in  the  design  of  the  thrust  vector  control 
law  (Chapter  III)  and  the  evaluation  of  the  controller  (Chapters  IV 
and  V ) .  The  equations  are . 


Lateral  Acceleration 


V  -  w0  p  +  u0  r  =  Y^fJ  +  Yg3  +  Ypp  +  Yrr  +  g  cos  0O  <j> 


+  Y.  6  +  Y.  6 

o  a  6  v 
a  v 


(1) 


Roll  Acceleration 

y  rr  • 

*  "  i-  =  L-0  +  Lg6  +  L  p  +  Lrr  +  6y 

a  v 


(2) 


Yaw  Acceleration 
Tv^  • 

r  -  1ST  j  °  V  +  V  +  V  +  V  +  "s  Sa  +  B6  {v 

a  v 


Kinematic 


<f>  =  p  +  r  tan  0O 


(3) 


(1*) 


The  quantities  v,  p,  r,<j,  and  3  represent  the  aircraft's  per¬ 
turbation  side  velocity,  roll  rate  ,  yaw  rate,  bank  angle  and  side¬ 
slip  angle  respectively.  The  subscript  "0"  represents  equilibrium 

values  and  the  quantities  Y„ ,  Y  ,  L  ,  etc.  are  the  dimensional 

p  P  P 

stability  and  control  derivatives  (Ref  2:6. U6 ) .  The  term  6  is  the 

a 

aircraft  perturbation  aileron  deflection  while  5^  represents  thrust 
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o 


vector  deflection.  The  <S^  terms  are  replaced  with  <5^  (rudder  deflec¬ 
tion)  in  the  tail-on  cases.  Appendix  3  gives  the  sicn  convention 

for  5  and  derives  the  TVC  dimensional  derivatives, 
v 

It  is  generally  assumed  that  the  8  derivatives  are  negligible 
(Ref  2: If. 82).  Then,  noting  that: 

V  =  vt  6  . 


and: 

V=vt8 


where  V  is  the  aircrafts  total  velocity. 
"0 

Equation  1  can  be  rewritten  as: 


.  Y  (w  +  Y  )  p  (Y  -  u  )  r  g  cos  0  <f» 

r - B  8  +  _ 2 _ £  +  _ r - o_  +  - a. 

p  "  V  V  V  V 

t  t  t  t 

Y  Y 

8  8  8  8 

,  a  a  v  v 


(5) 


Simultaneous  solution  of  Equations  (2)  and  (3)  for  p  and  r  and 
introduction  of  primed  dimensional  derivatives  into  these  equations, 
gives : 


P  =  Lo&  +  L"P  +  Kr  +  K*  K  +  Lx 

U  M  *  W  U  W  I 

r  V 

r  =  IO  +  K"p  +  N'r  +  Mi  6  +  Ni  6 

8  p1  r  6  a  6  v 

r  v 

The  primed  derivatives  are  defined  as: 


(6) 

(7) 


Lr 

1 


L. 

1 


1  - 


Ixz 

Izz 


Ixz 


2 


Ixx  Izz 


N. 

1 


8 


+  L. 

Izz  i 

-  Ixz2 
Ixx  Izr. 


where  the  subscript,  "i",  represents  B,  p,  r,  6  ,  6^. 

Finally,  the  lateral  acceleration  can  be  expressed  as: 


a  =  YaB  +  Y  p  +  Y,  6  +  Y  6 
y  3  p*  6&  a  6y  v 


(8) 
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III.  Thrust  Vector  Control  Lav  Development 


This  chapter  presents  the  control  lavs  as  developed  for  the  TVC 
yaw  stability  augmentation  system  (SAG).  The  requirements  for  the 
aircraft  dutch  roll  mode,  roll  mode  and  spiral  mode  are  presented. 

The  block  diagram  of  the  TVC  yaw  SAS  is  shown  and  compared  with  the 
•tail-on  yaw  SAS.  The  root  locus  plots  for  the  tail-off  aircraft  and 
yaw  SAS  are  shown  and  then  the  lateral  modes  for  the  closed  loop 
system  are  compared  with  the  requirements  and  the  tail-on  aircraft 
modes . 

Computer  Programs 

Existing  computer  programs  and  library  subprograms  available 
•through  the  Aeronautical  System  Division's  (ASD)  computer  system 
■were  utilized  in  the  TVC  control  law  development.  In  order  to  use 
•the  root  locus  analysis  method,  the  aircraft  lateral-directional 
transfer  functions  needed  to  be  obtained.  A  program  written  by 
J.  Li.  Griffin  (Ref  3)  accepts  the  lateral  non-dimensional  stability 
derfvatrves  as  input  and  computes  the  lateral-directional  transfer 
functions  in  both  polynomial  form  and  factored  form.  The  transfer 
functions  of  the  variables,  3,  p,  r  and  a  with  aileron  and  rudder 

y 

(or  thrust  -vector)  input  are  provided  as  output.  It  also  provides 
the  dutch  roll  frequency  and  damping  ratio,  roll  time  constant  and 
spiral  time  constant. 

Another  very  useful  computer  program,  Total,  (Ref  ^)  in  the  AFIT 
library  of  the  ASD  computer  system  was  used  for  the  root  locus  analysis. 
The  program  was  accessed  interactively.  Transfer  functions  were 
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provided  as  input  to  obtain  the  root  locus  dicital  printout  and  plots 
used  for  analysis  and  design. 

Requirements 

The  lateral-directional  flying  qualities  requirements  for  thr 
dutch  roll  mode,  roll  mode  and  spiral  stability  were  taken  from  the 
Air  Force  Flight  Dynamics  Laboratory  (AFFDL)  flying  qualities  speci¬ 
fication  (Ref  5:182,206,220).  The  requirements  are  specified  for  a 
class  IV  (highly  maneuverable)  aircraft,  category  A  (flight  phase  that 
requires  rapid  maneuvering  or  precise  flight  control),  level  1  (flying 
qualities  clearly  adequate  for  mission  accomplishment)  flying  qualities. 

The  three  lateral-directional  modes  (Ref  2:6.53)  can  be  obtained 
from  the  roots  of  the  characteristic  equation  for  the  aircraft.  The 
dutch  roll  mode  is  represented  by  a  complex  pair  of  roots  in  the 
s-plane  so  the  flying  qualities  requirements  are  specified  in  terms 
of  frequency  and  danping  ratio.  The  roll  mode  normally  appears  in  the 
characteristic  equation  as  a  root  on  the  s-plane  real  axis.  Its  re¬ 
quirement  is  specified  in  terms  of  the  roll  mode  time  constant  which 
implies  a  response  requirement  to  a  step  (aileron)  input.  The  spiral 
mode  is  often  divergent  so  the  requirement  is  specified  as  a  time  to 
double  amplitude.  The  requirements  for  these  three  modes  are  as 
follows : 

Dutch  Roll  Mode  -  w  >  ,  , 

n  -  1  rad/ sec 

C  -  0.19 

^Wn  -  0.35  rad/sec 
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1  per  sec 


Boll  Mode  -  1/t  > 

r  - 

Spiral  Stability-Minimum  time  to  double  amplitude  =  12  sec 
or  1/t  >  -  0.057  per  sec 

Coupled  Roll/Spiral  is  unacceptable. 

System  Analysi s 

By  looking  at  the  non-dimensional  stability  derivatives 

(Appendix  2)  it  can  be  seen  that  for  each  flight  condition  the 

stability  derivatives  are  different  for  the  tail-on  and  tail-off 

configurations.  The  stability  derivative  of  particular  interest 

here  is  C  since  this  derivative  must  be  positive  for  directional 

ng 

stability  (Ref  6:112-113).  Therefore,  it  can  be  seen  directly  from 
Appendix  2  that  removal  of  the  vertical  stabilizer  causes  the  air¬ 
craft  to  be  directionally  unstable  (C  is  negative).  The  effect 
of  the  removal  of  the  vertical  stabilizer  can  be  seen  from  the  large 

changes  in  C  ,  C  ,  C  and  C  .  In  order  to  restore  this  lost 
y„  y  n  n 

0  p  p  r 

directional  stability  a  feedback  control  system  to  a  thrust  vectoring 
device  is  to  be  used.  Figure  1  is  a  block  diagram  of  the  control 
system  proposed  to  correct  this  lack  of  directional  stability. 
Sideslip  (g)  feedback  is  used  in  the  inner  loop  with  yaw  rate  (r) 
feedback  in  the  outer  loop.  The  g  feedback  basically  increases  the 
dutch  roll  frequency  and  the  r  feedback  provides  the  dutch  roll  damp¬ 
ing. 

The  thrust  vector  actuator  was  selected  with  a  0.1  second  time 
constant.  This  was  believed  to  be  representative  of  an  actuator  for 
this  type  of  application..  The  washout  circuit  time  constant  was 
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Basic  Aircraft 


Washout 


selected  to  be  consistent  with  the  tail-on  washout  time  constant  of 


•1  second.  The  effect  of  washout  time  constant  will  be  discussed  later 
in  this  chapter.  A  simplified  block  diagram  of  the  tail-on  yaw  SAS 
is  shown  in  figure  2.  The  configuration  of  this  yaw  SAS  is  basically 
the  same  as  the  tail-off  yaw  SAS  except  that  the  lateral  acceleration 
(ar)  feedback  was  used  in  the  inner  loop.  The  rudder  actuator  is 
modelled  with  a  0.05  second  time  constant  so  it  is  a  faster  actuator 
than  the  one  selected  for  the  TVC  configuration. 

A  root  locus  of  the  inner  loop  (0  feedback)  of  the  tail-off  yaw 
SAS  at  Mach  0.8  is  shown  in  figure  3.  Of  particular  interest  here 
is  the  coupled  roll-spiral  mode  in  the  right-half  s-plane.  This 
condition  is  unacceptable  from  the  flying  qualities  requirements 
standpoint.  The  dutch  roll  poles  are  located  at  -0.60  +  0 . h 5 j  *  By 
increasing  the  feedback  gain  the  roll-spiral  roots  can  be  moved  back 
to  the  real  axis.  As  this  is  done  the  dutch  roll  mode  goes  unstable 
but  the  frequency  is  increased  above  the  minimum  requirement  of 
1  rad/sec.  Closure  of  the  yaw  SAS  outer  loop  stabilizes  the  dutch 
roll  as  shown  in  figure  4.  The  outer  loop  gain  is  selected  to 
maximize  the  product  of  the  dutch  roll  frequency  and  damping  ratio 
This  will  minimize  the  time  to  damp  the  transient  response. 

As  seen  in  figure  4  the  roll  and  spiral  modes  are  now  uncoupled 
(as  they  should).  The  gains  selected  were: 

Inner  Loop  =  16.22 

Outer  Loop  =  6.35 


Figure  3.  Yaw  SAS  '  Inner  Loop  Root  Locus  Fit.  Cond. 


Figure  4.  Yaw  SAS  Outei’  Loop  Root  Locus  Flt.Cond. 


The  closed  loop  transfer  function  is: 

r  _  -  M6  (s  +  .292) (s  +  .205  t  .931.1 ) (s  +  l) _ _____ 

%  "  (s  ♦  .M9)(.  -  .0U5)(s  ♦  .893  ±  2.06j)(s  ♦  U.ODU  ♦  5.5) 

This  can  be  compared  with  the  closed  loop  transfer  function  for 
the  Mach  0.8  tail-on  aircraft  (and  yaw  SAS): 

r  _  -  20. UU  (s  +  .268  ±  .7.1  )(s  +  .33^)(s  +  l) _ _ _ 

6  rp  (s  -  .0l6)(s  +  .HU6  ±  1.2,3 ) (s  +  .517)(s  +  2.83) (s  +  ik.3) 

Table  2  shows  the  dutch  roll  frequency  and  damping,  roll  time 
constant  and  the  spiral  time  constant  for  the  Mach  0.8  tail-off  case 
as  compared  with  the  flying  qualities  specification  requirements 
and  with  the  tail-on  case.  This  table  indicates  that  the  tail-off 
aircraft  exceeds  the  dutch  roll  specification  requirement.  The 
tail-off  configuration  has  a  higher  frequency  and  damping  than  the 
tail-on  configuration  which  suggests  a  better  transient  response 
than  the  tail  on.  The  spiral  divergence  criteria  is  met  for  the 
tail-off  aircraft  but  the  roll  time  constant  is  too  high.  Ho  attempt 
will  be  made  here  to  correct  this  problem  but  only  to  suggest  that 
roll  augmentation  would  correct  the  problem.  Another  way  to  correct 
the  problem  would  be  to  increase  the  washout  time  constant  sufficiently 
but  this  will  be  discussed  later  in  this  chapter. 


Table  2 


Mach  0.8  Lateral  Flying  Qualities  of  Interest 

Dutch  Roll  Roll 

„  „  1/ 

w  r,  sa>  t 

n  nr 

(rad/sec)  (rad/sec)  (l/sec) 

Requirement  1  (min)  0.29  (min)  0.35  (mini  1  (mini  -0^05?"  (min) 

Tail-off  2.2h  0.1*  0.89  0.1*7  -0.0l*5 

Tail-on  1.28  0.35  0.1*5  2.86  -0.0l6 

The  Mach  2  case  was  bandied  in  much  the  same  way  as  the  Mach  0.8 

case  but  new  feedback  gains  had  to  be  selected.  The  root  locus  of  the 

inner  loop  of  the  yaw  SAS  is  shown  in  figure  5*  This  time  the  tail-off 

aircraft  does  not  have  a  coupled  roll-spiral  or  a  dutch  roll  mode  but 

the  roots  are  all  on  the  real  axis.  By  increasing  the  feedback  gain 

•the  poles  at  2-9  and  0.022  form  the  dutch  roll  mode-  The  pole  at 

-3-15  couples  for  a  small  range  of  gain  with  the  pole  at  -1 . 1*3  but 

then  they  return  to  the  real  axis  to  form  a  roil  mode  and  spiral  mode. 

The  dutch  roll  is  unstable  for  all  values  of  gain  but  the  yaw  rate 

feedback  will  stabilize  this  mode.  An  inner  loop  gain  vas  selected 

which  produced  a  good  dutch  roll  root  locus  (figure  6).  The  gains 

selected  were: 

K  -  6.1*3 
8 

Kr 


Spi  ral 


(l/sec) 
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Figure  5-  Yaw  SAS  Inner  Loop  Root  Locus  Fit-  Cond. 


The  transfer  function  for  the  closed  loop  system  (figure  l)  is: 

r  _  -  .32.8  (s  +  . 862 ) ( s  +  .335  ±  .6.1  )(s  +  1) _ 

6r?  (s  +  .96)(s  -  .0108 )(s  +  1.24  ±  2.76j)(s  +  4.03)(s  +  5-2) 

For  comparison  the  transfer  function  for  the  tail-on  aircraft 
vith  its  yaw  SAS  (figure  2)  is: 

r  _  -  42  (s  +  .65  t  1 ,05.i )  (s  +  .65)  (s  +  l) _ 

* 

r  (s  +  . 009 ) ( s  +  2.12  ±  6.17,5 )(s  +  -94)(s  +  2.77)(s  +  7-37) 

P 

Table  3  shows  the  dutch  roll  frequency  and  damping,  roll  time 
constant  and  the  spiral  time  constant  for  the  Mach  2  tail-off  case  as 
compared  vith  the  flying  qualities  specification  requirements  and  the 
Mach  2  tail-on  case.  While  the  tail-off  aircraft  exceeds  the  spec¬ 
ification  requirement  ,  the  tail-on  aircraft  has  a  much  higher  fre¬ 
quency  and  parameter  which  indicates  that  it  will  have  a  shorter 
transient  response  than  the  tail-off  aircraft.  The  roll  and  spiral 
time  constants  satisfy  the  specification  requirements. 


Table  3 

Mach  2  Lateral  Flying  Qualities  of  Interest 

Dutch  Roll  Roll  Spiral 


w 

n 

(rad/sec) 

C 

n 

(rad/sec) 

H 

r 

( 1/sec) 

r 

(l/sec) 

Requirement 

1  (min) 

0.19  (min) 

0.35  (min) 

1  (min) 

-0.057  (min) 

Tail-off 

3.0 

o.ia 

1.24 

4.0 

-0.0108 

Tail-on 

6.5 

0.32 

2.12 

2.78 

0.009 
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Effects  of  Washout  Time  Constant 


As  indicated  previously  the  washout  circuit  for  the  tail-off 
design  was  selected  to  be  consistent  with  that  in  the  tail-on  yaw 
SAP-.  The  effect  on  the  dutch  roll  of  changing  from  the  1  second 
time  constant  was  examined  for  flight  condition  1.  The  result  is 
that  as  the  time  constant  is  decreased  the  root  locus  of  the  dutch 
roll  collapses  toward  the  imaginary  axis.  Figure  7  shows  the  effect 
of  0.83  second  washout  time  constant.  If  the  washout  time  constant 
is  increased  the  result  is  that  the  dutch  roll  root  locus  expands 
away  from  the  imaginary  axis  producing  much  better  dutch  roll 
characteristics.  Figures  8  and  9  show  the  effect  of  a  1.25  second  and 
a  2  second  time  constant  respectively. 

If  a  washout  time  constant  on  the  order  of  2  seconds  could  be 
tolerated,  the  washout  pole  would  be  to  the  right  of  the  roll  pole 
(figure  9)  and  the  branch  associated  with  the  roll  mode  would  move 
to  the  left.  There  would  then  be  no  problem  in  satisfying  the  flying 
qualities  roll  time  constant  requirement. 
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igure  7.  Effect  of  a  .83  sec  Vasftout  Tiir.e  Constant  Fit.  Ccnd. 


Figure  9-  Effect  at  a  2  sec.  V.'ashout  Time  Constant  Fit.  Cond. 


at.  a  subsonic  cruise  condition.  However,  it  takes  the  tail-on  aircraft 
longer  to  damp  the  disturbance.  This  would  be  expected  as  the  dutch 
roll  frequency  and  damping  are  less  than  that  of  the  tail-off  air¬ 
craft  at  this  flight  condition.  Figure  19  shows  plots  of  p  and  $  for 
the  tail~on  aircraft.  The  bank  angle  seeks  a  steady  state  value  of 
about  9  degrees  whereas  the  tail-off  aircraft  continued  to  roll,  but, 
at  a  slow  rate.  After  10  seconds  the  tail-off  aircraft  only  rolled  3 
3  degrees.  This  reflects  the  long  roll  mode  time  constant. 

The  Mach  2  tail-off  rudder  pulse  responses  are  shown  in  figures 
20-25.  Figure  20  shows  the  pulse  rudder  input,  6^,  and  6^  plots. 

The  variables  a  3  and  r  plots  are  shown  in  figure  21.  Again,  the  yaw 

y  j 

SAS  is  functioning  properly  although  this  time  the  magnitude  of  a  ,  3 

y 

and  r  are  larger  since  the  thrust  is  larger  than  the  Mach  0.8  case. 

The  response  damps  out  quickly  (less  than  k  sec.)  which  would  be 
expected  from  the  high  dutch  roll  frequency  and  damping  ratio. 

The  plots  of  p  and  4>  are  shown  in  figure  22.  In  this  case  the  bank 
angle  reached  a  steady  value. 

Figure  23  shows  the  Mach  2  tail -on  a  ,  3  and  r  pulse  rudder 

y 

plots.  The  magnitudes  of  these  variables  are  not  as  high  as  in  the 
tail-off  case  and  the  system  damps  out  faster  (3  sec.).  This  would  be 
expected  as  the  product  of  the  dutch  roll  frequency  and  damping  ratio 
is  higher.  The  rudder  position  and  rate  plots  are  shown  on  figure  2k 
and  reflect  the  high  dutch  roll  frequency.  The  higher  a  in  the 

y 

tail-off  case  reflects  the  high  thrust  at  this  flight  condition. 

The  plots  of  <J>  and  p  are  shown  in  figure  25.  For  x  'light  case 
the  bank  angle  reached  a  lower  steady  state  va  .n  the  tail-off 
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(Ref  5:M8)  were  used  to  calculate  the  appropriate  value  for  v  .  In 

y 

order  to  get  a  P  gust  term  (P  )  the  following  relation  is  used: 

6 


f*  =  -  ^ 


The  discrete  P^  is  then  input  through  the  P  term  in  the  equations  of 
motion. 


Random  Turbulence 

To  complement  the  evaluation  of  the  aircraft  and  control  system 
to  atmospheric  turbulence,  the  effect  of  continuous  random  turbulence 
needs  to  be  considered.  The  reason  is  that  if  a  measure  of  the  velo¬ 
city  of  gusts  on  a  given  flight  path  were  made,  the  magnitude  of  the 
gusts  would  be  of  a  random  nature.  Therefore,  continuous  random 
turbulence  is  described  statistically.  Again  only  side  gusts  will 
be  considered  since  the  study  is  limited  to  the  lateral-directional 
modes . 

In  straight,  steady  level  flight  through  a  turbulent  atmosphere 
it  is  important  to  minimize  the  unwanted  motion  (error)  caused  by 
the  turbulence.  In  order  to  statistically  measure  how  well  the  SAS 
is  performing,  mean  square  error  is  utilized  since  the  average  value 
of  the  gusts  is  assumed  to  be  zero.  The  power  spectral  density,  $ , 
is  utilized  to  obtain  the  mean  square  error  for  an  aircraft  and 
control  system.  The  random  turbulence  model  is  represented  in 
figure  10.  The  input  used  to  evaluate  the  performance  of  the  system 
is  white  noise  which  has  a  constant  power  spectral  density  (=  l). 

There  are  several  textbooks  which  give  a  treatise  on  random  atmospheric 
turbulence.  Blakelock  is  one  (Ref  6:  Chapter  9). 
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The  Dryden  form  of  the  power  spectral  density  for  the  side  gust 
is  utilized  and  is  defined  as : 


1  +  3  ( 

$>  (u)  =  a  Lv  _ _t _ 

S  Vt  1  +  [,  Lv  ^21  2 

.  t 


where  Lv  and  are  parameters  specified  in  reference  5  (p  M7-M8), 
The  transfer  functions  for  the  shaping  filter  (T  (j  a)  )  can  be 

vs 

obtained  by  performing  spectral  factorization  of  <I>  (to).  That  is: 

Vg 


v  (oj)  =  T  (jw)  T  (-jw) 
6  g  g 


1+/3  ju) 

\  (j»)  =  %  HE  _ t 

^  (it  0W)2 

t 


Then  by  replacing  jw  with  s  and  dividing  by  V,  yields. 


T  (s)  =  a  I Lv 


1+  ]3  £  s 


(1  +  I1  s)2 
t 


■which  is  utilized  in  the  continuous  random  turbulence  model. -The 
values  (thunderstorm  levels)  for  Lv  and  used  are: 

Lv  =  1750  ft. 


Cy.  =  21  ft/sec 
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Computer  Models /Programs 

In  order  to  obtain  time  response  plots  for  the  pulse  rudder  pedal 
input  and  the  discrete  atmospheric  turbulence  input,  an  integration 
and  plot  computer  program  was  created.  The  integration  was  done  using 
a  library  Runge-Kutta  integration  subroutine  (ODE)  on  the  CC6600 
library  of  the  ASD  computer  system.  The  aircraft  and  control  system 
modeled  in  state  variable  form  were  provided  as  input  to  this  program. 
The  form  of  the  program  input  is: 

t  =  A  X  +  G  S  +  B6 


where  X  is  the  aircraft  state  vector,  6  is  pilot  rudder  pedal  inputs, 

r  * 

P 

and  Bg  is  gust  inputs. 

Figure  11  shows  the  block  diagram  of  the  state  space  model  used 
for  the  tail-off  configurations.  The  only  differences  in  the  model 
for  the  tail-on  aircraft  are  the  a  feedback  and  the  different 

y 

actuator  dynamics.  The  state  equations  in  matrix  form  are  shown  in 
figure  12.  The  gust  term  (B^)  enters  each  of  the  state  equations 
through  the  B  term.  For  the  pilots  rudder  pedal  input  the  G  matrix 


is  a  null  matrix.  For  the  discrete  gust  input  the  B  matrix  is  a  null 


matrix.  The  output  of  the  computer  program  is  time  response  plots 

« 

(0-10  second)  for  $  +  f$  ,  p,  r,  i{> ,  6  ,  6  ,  a  ,  and  ft  (or  6  )  and 

S  v  v  y  g  r^ 

a  listing  at  0.05  second  increments  for  these  variables. 


An  existing  computer  program  developed  by  Capt  J.  T.  Silverthorn 
(AFIT  Department  of  Aeronautics  and  Astronautics),  which  performs  a 
covariance  analysis  (Ref  7:89)  was  utilized  to  evaluate  the  response 
of  the  aircraft  to  random  atmospheric  turbulence.  The  inputs  to  the 
program  are  the  aircraft  and  control  system  state  equations  augmented 
with  the  state  equations  for  the  shaping  filter  to  model  the  atmo¬ 
spheric  power  spectral  density  function.  The  power  spectral  density 
function  adds  two  states  to  the  basic  aircraft  and  control  system 
A  matrix  (figure  12).  A  listing  of  this  program  is  shown  in  Appendix  5. 

Response  To  A  Rudder  Pedal  Input 

In  order  to  determine  the  effectiveness  of  thrust  vectoring  to 
generate  a  yawing  moment,  a  comparison  with  the  tail-on  rudder 
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Figure  12.  Integration  Program  Input  Matrix 


effectiveness  vas  made.  .This  was  expressed  in  terms  of  yawing 

moment /degree  of  thrust  vector  deflection  (izs  IL  )  compared  with 

v 

the  yawing  moment/degree  of  rudder  deflection  (izz  N.  ).  This  is 

6 

r 

plotted  in  figure  13  for  different  Mach  numbers  at  3Q000  ft,  altitude. 

What  is  indicated  is  that  for  the  subsonic  Mach  numbers  the  rudder 

generates  a  higher  control  yawing  moment  than  does  thrust  vectoring. 

Note  the  low  thrust  condition  at  Mach  0.8.  But  when  the  thrust  is 

increased  sufficiently  in  the  supersonic  regime,  the  thrust  vector 

generates  a  higher  control  yawing  moment. 

A  command  rudder  (6  )  pxilse  of  3  degrees  for  1  second  duration 

rp 

vas  applied  to  the  system  to  observe  the  response.  In  figure  1^  is  a 
0-10  sec.  plot  of  the  input  and  plots  of  the  position  of  the  thrust 
vector  (6v)  and  rate  which  the  thrust  vector  is  moving  (6^)  for  the 
Mach  0.8  tail-off  case.  It  can  be  seen  from  the  6  plot  that  the 
yaw  SAS  is  resisting  the  input  as  it  should  do  initially.  Figure  15 
shows  plots  of  e  ,  P  and  r  where  the  dutch  roll  mode  can  be  seen. 

y 

Notice  that  a  has  the  same  shape  as  6  which  indicates  that  the 
y  v 

lateral  acceleration  is  primarily  the  result  of  thrust  vectoring.  The 
amplitudes  of  a  and  B  are  small  and  the  time  to  damp  out  is  less  than 
5  seconds.  Figure  l6  shows  plots  of  p  and  <j>. 

The  same  rudder  pulse  was  applied  to  the  tail-on  aircraft  to 
compare  its  response  to  that  of  the  tail-off  aircraft.  Figure  17 
shows  plots  of  the  rudder  position  and  rate.  Figure  18  shows  plots 
of  a  ,  B  and  r  which  indicate  larger  amplitudes  than  for  the  tail-off. 

y 

This  reflects  what  was  shown  in  figure  13  that  the  rudder  control 
effectiveness  is  greater  than  the  thrust  vector  control  effectiveness 
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at.  n  subsonic  cruise  condition.  However,  it  takes  the  tail-on  aircraft 
longer  to  damp  the  disturbance.  This  would  be  expected  as  the  dutch 
roll  frequency  and  damping  are  less  than  that  of  the  tail-off  air¬ 
craft  at  this  flight  condition.  Figure  19  shows  plots  of  p  and  <j>  for 
the  tail~on  aircraft.  The  bank  angle  seeks  a  steady  state  value  of 
about  9  degrees  whereas  the  tail-off  aircraft  continued  to  roll,  but, 
at  a  slow  rate.  After  10  seconds  the  tail-off  aircraft  only  rolled  3 
3  degrees.  This  reflects  the  long  roll  mode  time  constant. 

The  Mach  2  tail-off  rudder  pulse  responses  are  shown  in  figures 
20-25.  Figure  20  shows  the  pulse  rudder  input,  6^,  and  6^  plots. 

The  variables  a  B  and  r  plots  are  shown  in  figure  21.  Again,  the  yaw 

y  > 

SAS  is  functioning  properly  although  this  time  the  magnitude  of  a  ,  B 

y 

and  r  are  larger  since  the  thrust  is  larger  than  the  Mach  0.8  case. 

The  response  damps  out  quickly  (less  than  sec.)  which  would  be 
expected  from  the  high  dutch  roll  frequency  and  damping  ratio. 

The  plots  of  p  and  $  are  shown  in  figure  22.  In  this  case  the  bank 
angle  reached  a  steady  value. 

Figure  23  shows  the  Mach  ?  tail -on  a  .  B  and  r  pulse  rudder 

y 

plots.  The  magnitudes  of  these  variables  are  not  as  high  as  in  the 
tail-off  case  and  the  system  damps  out  faster  (3  sec.).  This  would  be 
expected  as  the  product  of  the  dutch  roll  frequency  and  damping  ratio 
is  higher.  The  rudder  position  and  rate  plots  are  shown  on  figure  2h 
and  reflect  the  high  dutch  roll  frequency.  The  higher  a  in  the 

y 

tail-off  case  reflects  the  high  thrust  at  this  flight  condition. 

The  plots  of  $  and  p  are  shown  in  figure  25.  For  u  'light  case 
the  bank  angle  reached  a  lower  steady  state  va  .n  the  tail-off 


aircraft  due  to  the  shorter  time  to  damp  the  dutch  roll  mode. 

Response  To  Discrete  Atmospheric  Turbulence 

In  order  to  properly  evaluate  the  response  of  the  aircraft  and 
SAS  to  discrete  (l  -  COS)  atmospheric  turbulence,  the  gust  frequency 
which  excites  the  aircraft  the  most  (resonance)  must  be  found. 

A  log  magnitude  frequency  response  (Bode  plot)  of  the  system  to  a 
sinusoidal  input  could  be  useful  in  finding  the  resonant  frequency. 
Noting  the  half  angle  formula: 


sin  f<» 

>1  + 

1  -  cos  <*> 

r 

c 

r 

2 

Where:  u>  =  frequency 

a  relationship  of  a  sinusoidal  input  to  a  1  -  COS  input  can  be  seen. 
Basically,  th  resonant  frequency  of  a  sinusoidal  input  is  that  of 
the  1  -  COS  input.  Frequency  response  plots  of  g/  g  and  a  /g 

s  y 

are  shown  ir.  figures  2 6  and  27  for  the  Mach  0.8  tail-off  it 
case.  The  resonant  frequency  is  shown  to  be  about  2  rad/seo.  It 
might  be  expected  that  the  resonant  frequency  for  a  1  -  COS  discrete 
input  would  be  about  Ij  rad/sec. 

A  frequency  sensitivity  study  was  done  using  the  integration 
program  described  earlier  in  this  chapter.  The  results  indicated 
that  a  gust  frequency  of  3  rad/sec  produced  the  largest  values  of 
g  and  a  (Table  ll).  A  gust  frequency  of  ^  rad/sec  excited  the  TVC, 

y 

with  the  highest  <$v  excursion.  For  the  tail-on  aircraft  a  gust 
frequency  of  2  rad/sec  produced  the  largest  values  of  g  and  a  . 

y 

From  Table  it  is  seen  that  as  the  gust  frequency  increases  the 
peak  thrust  vector  excursion  increases  until  a  b  rad/sec  gust  is 
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Discrete  Gust  Frequency  Sensitivity  Study 


reached  "but  the  rate  continues  to  climb  with  gust  frequency.  The 
same  is  true  for  the  tail-on  rudder  excursion  and  rate  but  the 
magnitudes  are  much  smaller. 

The  frequency  response  plots  for  the  Mach  2  case  are  shown  in 
Figures  28  and  29.  The  resonant  frequency  to  a  sinusoidal  input 
is  2.5  rad/sec.  A  frequency  sensitivity  study  using  the  discrete 
gust  input  was  also  done  for  the  Mach  2  flight  condition,  Table  5 
shows  the  results.  For  this  case  higher  frequencies  excited  the 
system  since  the  dutch  roll  frequencies  for  both  tail-on  and 
tail-off  were  higher  than  in  the  Mach  0.8  case.  For  the  Mach  2 
tail-off  case  a  gust  frequency  of  5  rad/sec  produced  the  largest 
values  of  B  and  a^,.  The  largest  value  for  6^  was  also  produced  at 
this  frequency  but  <5^  increased  with  the  gust  frequency.  The  tail- 
on  aircraft  has  a  higher  dutch  roll  frequency  (6.53  rad/sec)  than 
the  tail-off  aircraft  so  the  largest  values  of  3  and  a  occurred 

y 

at  a  higher  gust  frequency  of  6  rad/sec.  For  this  flight  case 
the  values  for  the  rudder  excursion  and  rates  were  comparable  to  the 
thrust  vector  excursion  and  rates.  This  would  be  expected  as  the 
yawing  moment  due  to  thrust  vectoring  (figure  13)  is  much  higher  at 
the  Mach  2  flight  condition. 

Time  response  plots  were  run  for  the  Mach  0.8  flight  condition 

at  a  gust  frequency  of  ^  rad/sec.  This  frequency  was  chosen  as  it 

most  excited  the  TVC,  Time  response  plots  for  the  tail-off  aircraft 

are  shown  in  figure  30-32  and  the  tail-on  plots  are  shown  in  figures 

33-35*  From  figure  30  the  3  input  can  be  seen  and  also  the  6  and 

g  v 

6^  time  responses.  Figure  31  shows  the  a  ,  3  and  r  time  responses. 
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In  comparing  a  and  g  parameters  to  those  of  the  tail-on  aircraft 

y 

responses  in  figure  34  it  is  seen  that  the  tail-off  case  produces 
larger  values  of  (3  and  a  hut  damps  out  faster.  Since  the  tail-off 

y 

aircraft  has  a  larger  dxitch  roll  frequency  and  damping  than  the 
tail-on  aircraft  the  faster  damping  would  he  expected.  Looking  at 
the  and  time  responses  one  realizes  the  price  which  must  he  paid 
for  this  stability.  peaks  at  38  degrees  and  6^  peaks  at  119 

deg/sec  (as  indicated  in  table  4).  In  comparing  the  <j>  and  p 
parameters  in  figures  32  (tail-off)  and  35  (tail-on)  the  same  type 
of  transient  behavior  is  disp3.ayed  as  the  other  aerodynamic  parameters 
except  for  a  slightly  larger  hank  angle  in  the  tail-on  case  than  in 
the  tail  off  case. 

For  the  Mach  2  flight  condition  a  discrete  gust  frequency  of 
5  rad/sec  was  used  since  this  excited  the  tail-off  aircraft  the 
most.  This  was  also  anticipated  from  the  frequency  response  study. 
Time  response  plots  for  the  tail-off  aircraft  are  shown  in  figures 

36-38.  The  input  can  he  seen  in  figure  36  along  with  <$v  and 

•  » 

<5^  and  it  is  immediately  noticed  that  the  magnitudes  of  6^  and  6^ 

are  much  lower  and  more  reasonable  than  the  Mach  0.8  flight  condition. 

From  figure  37  the  transient  response  is  seen  to  damp  much  faster 

than  in  the  Mach  0.8  case  due  to  the  higher  dutch  roll  frequency  and 

damping.  The  high  dutch  roll  frequency  (6.5  rad/sec)  of  the  tail-on 

aircraft  can  he  seen  in  figures  39  and  40  as  the  transient  response 

is  of  high  frequency  and  damps  out  quicker  than  the  tail-off 

aircraft.  The  tail-off  aircraft  provides  a  softer  response  than  the 

tail-on  aircraft  as  the  lateral  acceleration  is  lower  although  P  is 

slightly  larger.  Figure  hi  shows  the  time  responses  of  and  p. 
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Response  To  Continuous  Random  Atmospheric  Turbulence 


The  response  of  the  aircraft  and  control  system  to  continuous 
random  turbulence  needs  to  be  seen  in  order  to  complement  the 
discrete  turbulence  evaluation.  The  random  turbulence  model  and 
computer  program  explained  previously  were  utilized  for  this  part 
of  the  study.  The  input  is  in  the  form  of  white  noise  and  the  output 
is  a  matrix  of  variances  or  mean  square  error  values  for  each  state 
variable  about  a  mean  of  zero.  The  results  are  presented  in  table  6 
as  the  square  root  of  the  variances  or  standard  deviation  of  the  error. 

From  table  6  it  is  seen  that  for  the  Mach  0.8  flight  condition 
the  standard  deviation  of  the  aircraft  output  variables  are  reasonably 
close  comparing  the  tail-off  and  tail-on  values.  Once  again  rel¬ 
atively  large  thrust  vector  excursions  and  rates  are  seen.  In  the 
Mach  2  flight  condition  the  tail-off  aircraft  actually  has  a  smaller 
standard  deviation  than  the  tail-on  aircraft  due  to  the  large  amount 
of  yawing  moment  (figure  13)  available.  The  tail-on  aircraft  exhibits 
a  larger  lateral  acceleration  at  Mach  2  as  it  did  for  the  discrete 
turbulence . 

The  results  shown  in  this  chapter  have  indicated  that  high  thrust 
vector  excursions  and  rates  are  commanded  by  the  SAS  in  a  low  thrust 
flight  condition  (flight  condition  l).  While  the  linear  assumptions 
made  in  the  beginning  become  invalid,  the  trends  should  remain  un¬ 
changed.  The  next  step  is  to  limit  the  thrust  vector  excursion  and 
rate  to  more  realistic  values.  The  next  chapter  will  show  these 
effects.  This  feature  also  allows  the  linear  assumptions  to  be  accept¬ 
able  so  that  the  small  perturbation  equations  can  still  be  used. 
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V.  Thrust  Vector  Control  Requirements 


From  the  last  chapter  it  was  seen  that  large  values  of  thrust 
vector  excursion  and  rate  were  commanded  "by  the  control  system  due 
to  atmospheric  turbulence  r.t  the  thunderstorm  levels.  Practically, 
there  are  limits  imposed  by  the  hardware.  In  addition  longitudinal 
effects  become  increasingly  important  when  the  TV  excursion  is  large. 
In  this  chapter,  the  results  of  a  study  of  the  effects  of  limiting 
the  thrust  vector  excursion  and  rate  is  shown  in  order  to  establish 
the  basic  requirements  to  provide  acceptable  stability.  By  limiting 
the  TV  excursion  the  small  perturbation  equations  of  motion  are 
again  utilized  with  the  longitudinal  effects  being  ignored. 

This  chapter  also  examines  the  sensitivity  of  the  closed  loop 
system  to  variations  in  thrust,  feedback  gains,  and  actuator  dynamics. 

The  effect  on  the  lateral-directional  modes  of  utilizing  a 
fast  and  a  slow  TVC  actuator  is  also  addressed. 

Thrust  Vector  Limit  Study 

The  thrust  vector  excursion  and  excursion  rate  was  limited  to 
see  the  effect  on  the  dutch  roll  frequency  and  damping.  The  discrete 
atmospheric  turbulence  model  was  again  utilized.  The  integration 
program  (Appendix  h)  was  modified  to  permit  limiting  of  both  the 
thrust  vector  deflection  and  rate  of  deflection.  As  the  limits  were 
decreased  the  point  where  the  system  did  not  meet  the  flying  qualities 
requirements  was  determined.  Further  decrease  in  the  limits  caused 
the  aircraft  to  become  unstable. 

As  seen  previously  the  maximum  thrust  vector  excursion  for  a 


discrete  side  gust  of  h  rad/sec  was  38  degrees  with  a  maximum  ex¬ 
cursion  rate  of  119  deg/sec  (Table  U).  A  limit  on  the  TV  excursion 
of  0.5  radians  (28.7  degrees)  with  an  excursion  rate  limit  of  1.2 
rad/sec  (68.8  deg/sec)  was  tried  first.  Then  subsequent  reductions 
in  the  rate  limit  (holding  the  TV  excursion  limit  constant)  were 
made  until  divergence  of  the  dutch  roll  mode  was  seen.  This  process 
of  decreasing  the  rate  while  holding  the  excursion  constant  was  con¬ 
tinued  for  lower  values  of  TV  excursion  do™  to  0.2  radians 
(11.5  degrees)  of  a<..f lection. 

In  Table  7A  a  tabulation  of  the  8  peaks  for  various  TV  excursion 
and  rate  limits  are  shown.  The  response  peaks  at  each  half  cycle 
of  the  dutch  roll  mode  are  sho™  and  their  time  of  occurance.  In 
Table  7B  is  sho™  the  dutch  roll  parameters  of  merit  for  each  limit 
case.  These  parameters  were  obtained  from  the  time  and  subsidence 
ratio  for  the  second  and  fourth  g  peaks  given  in  Table  7A.  It  can 
be  seen  from  the  table  that  decreasing  the  TV  rate  decreases  the 
dutch  roll  frequency  and  damping.  In  general,  if  the  TV  excursion 
is  decreased  with  the  rate  held  constant  the  dutch  roll  frequency 
decreases  but  the  damping  increases.  The  effects  of  holding  the 
thrust  vector  maximum  excursion  at  0.5  radians  and  decreasing  the 
maximum  rate  from  1.2  rad/sec,  to  1  rad/sec,  and  then  to  0.9  rad/sec 
are  sho™  in  figures  U2— 50 .  For  the  last  case  the  dutch  roll  mode 
is  divergent.  Seen  from  the  plots  is  that  as  the  TV  rate  is  decreased 
the  dutch  roll  mode  takes  longer  to  damp  out.  From  the  6  and  <5^ 
plots  it  is  seen  that  as  the  allowable  6  is  decreased  (slower  moving) 
the  thrust  vector  stays  at  the  deflection  limit  longer.  Eventually 
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(Fit.  Cond.  1,  U  rad/sec  Discrete  Gust) 


r 


o 


1.85  3.56  3.95  -i.su 


Table  7  B 


6  and  6  Limit  Study 
v  v 

(Fit.  Cond.  1,  4  rad/sec  Discrete  Gust) 

Dutch  Roll  Parameters 


6 

V 

(rad) 

6 

V 

(rad/sec) 

0) 

n 

(rad/sec) 

Co) 

n 

(rad/sec) 

Tail-on 

1.28 

0.35 

0.45 

Tail-off 

0.66 

2.08 

2.20 

0.35 

0.77 

0.5 

1.2 

2.17 

0.26 

0.56 

0.5 

1.0 

1.85 

0.085* 

0.16s 

0.4 

1.2 

2.09 

0.28 

0.59 

0.4 

1.0 

1.90 

0.13* 

0.25* 

0.3 

1.2 

1.90 

0.28 

0.53 

0.3 

1.0 

1.80 

0.22 

0.40 

0.25 

1.2 

1.88 

0.29 

0.55 

0.25 

1.0 

1.74 

0.23 

0.4o 

0.25 

0.9 

1.70 

0.20 

0.34* 

0.2 

1.2 

1.75 

0.28 

0.49 

0.2 

3  .1 

1.75 

0.26 

0.46 

*  Does  Not  Meet  Flying  Qualities  Specification 
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Figure  Mi.  Fit.  Cond.  1,  Tail-off 
Discrete  Atmospheric  Turbulence 

v  llax  =  .5  rad,  Max  =  1.2  rad/sec 
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Figure  ^8.  Fit.  Cond.  1,  Tail-off,  Discrete 
Atmospheric  Turbulence 
6  Max  =  .5  rad,  <$  Max  =  .9  rad/sec 


1,  Tail-off,  Discrete  Atrospheric  Turbulence 
•  5  rad,  <5^  Max  =  .9  rad/sec 


it  can  no  longer  stabilize  the  aircraft.  From  tables  7  A  and  7B  it 
can  also  be  seen  that  if  the  TV  excursion  is  reduced  the  rate  can  be 
reduced  also  (to  a  point)  and  still  maintain  acceptable  stability. 

To  further  substantiate  this,  time  response  plots  for  a  maximum  TV 
deflection  of  0.3  radians  and  a  maximum  rate  of  1  rad/sec  are  shown  in 
figures  51-53.  From  these  plots  it  is  seen  that  for  a  constant  TV 
excursion  rate,  the  dutch  roll  mode  damps  out  faster  than  for  the  case 
with  a  larger  excursion  (figures  ^5-^7). 

The  limit  study  can  be  summed  up  in  figure  5^  in  which  satis¬ 
factory  combinations  of  TV  deflection  and  TV  rate  of  deflection 
for  the  discrete  atmospheric  turbulence  environment  (side  gust)  is 
shown.  Areas  of  convergence,  marginal  stability  and  divergence  of 
the  dutch  roll  mode  are  indicated.  This  plot  shows  the  deflection 
and  rate  requirements  to  meet  the  flying  qualities  specification 
requirements  for  the  dutch  roll  mode. 

No  attempt  was  made  to  limit  the  thrust  vector  maximum  deflection 
and  rate  for  flight  condition  2  as  the  maximum  values  were  low 
(table  5).  It  is  believed  that  flight  condition  1  will  drive  the 
design  choice  (for  cruise). 

Parameter  Study 

A  study  of  the  effect  on  the  dutch  roll  frequency  and  damping 
of  varying  the  thrust,  feedback  gains  and  TVC  actuator  dynamics  by 
+  10$  was  done.  The  TV  maximum  deflection  was  limited  to  0.5  radians 
and  the  rate  was  limited  to  1.2  rad/sec.  Both  the  discrete  and  the 
continuous  random  turbulence  models  were  used  for  this  study. 


OELTft  V  -  DEO/SEC  OEl.TR  V  -  OEO  BETR  GUST 


(Rad/Sec) 
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Table  8  indicates  the  results  for  discrete  atmospheric  tur¬ 
bulence  at  flight  condition  1.  The  Maximum  values  of  0  and  a  are 

y 

tabulated  along  with  the  dutch  roll  mode  parameters.  What  is  in¬ 
dicated  is  that  the  system  is  most  sensitive  to  the  actuator  time 
constant  variations.  For  the  slower  actuator  the  dutch  roll  mode 
became  unacceptable  from  the  requirements  standpoint. 

The  results  of  the  parameter  study  for  flight  condition  2  are 

shown  in  table  9.  Values  for  6  and  6  were  included  since  the 

v  v 

magnitudes  of  these  quantities  were  not  limited.  It  can  be  seen 
from  table  9  that  the  system  is  less  sensitive  to  the  actuator 
variations  than  to  the  thrust  and  gain  changes.  What  is  indicated 
is  that  the  gains  need  to  be  closely  correlated  with  the  thrust 
and  the  actuator  needs  to  be  accurately  modelled.  Table  10  shows 
the  effects  of  the  parameter  variations  using  the  continuous  random 
atmospheric  turbulence  input.  It  substantiates  the  previous  results. 

TVC  Actuator  Model 

The  basic  TVC  actuator  model  used  as  a  baseline  in  this  study 
utilized  a  0.1  second  time  constant.  It  could  be  useful  from  a  hard¬ 
ware  standpoint  to  see  how  slow  an  actuator  could  be  used  and  still 
maintain  acceptable  flying  qualities.  This  would  be  of  interest  in 
the  selection  of  the  type  of  TVC  device  to  be  used.  Values  of  the 
actuator  time  constant  were  tried  at  0.2,  0.25  and  0.33  seconds. 

The  design  process  described  in  Chapter  III  was  repeated  to  develop 
new  feedback  gains  to  optimize  the  dutch  roll  dynamics.  It  was 
determined  that  with  an  actuator  time  constant  larger  than  0.25 
seconds  the  coupled  roll-spiral  mode  in  flight  condition  1  would 
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Table  8 


Parameter  Study,  Fit.  Cond.  2 
Discrete  Atmospheric  Turbulence 


Peak  Talues 

3  ay 

*  2 

(deg)  (ft/sec) 

Dutch  Boll 

a  l 

n 

(rad/sec) 

n 

(rad/sec) 

Tail-On 

2.25 

1.5 

1-28 

0.35 

0-.it5 

Tail-off 

Baseline 

3.98 

2-i t 

2-17 

0.26 

0-56 

-lO/i  Thrust 

3.76 

2.2 

1-96 

a.  32 

0,63 

-tlo#  Thrust 

it  .17 

2..  6 

2-2it 

Q.18* 

0-i<0 

9 

Actuator  — ^ 

it.  1(9 

2.5 

1.96 

D-D  8* 

0  .16* 

Actuator  ", 
s+11 

3.58 

2..  3 

2-03 

0-32 

D-65 

—10$  Gains 

3.80 

2,-lt 

2- 09 

0-28 

0  .59 

-+10$  Gains 

It.  lit 

2  .it 

2-17 

0.22 

0  -18 

■V  Max  -  .5  rad 

"v  Max  =  hh  1.2  rad/sec 

*  Does  Not  Meet  Flying  Qualities  SpeciTi cation 
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Parameter  Study,  Fit.  Cond. 


*uj 


+105?  Gains  i  9-3  1.82  10.58  *40.90  |  3.*4*4 


Table  10 


not  uncouple,  which  is  an  unacceptable  condition.  With  an  actuator 
time  constant  of  0.25  seconds  an  acceptable  dutch  roll  frequency  and 
damping  ratio  was  obtained  although  it  was  not  as  good  as  the  baseline 
actuator.  Table  11  is  a  tabulation  of  the  figures  of  merit  for  the 
lateral-directional  modes  for  the  baseline  actuator,  the  slower 
actuator  and  a  0.05  second  time  constant  (faster)  actuator.  The 
faster  actuator  was  included  to  indicate  its  advantage  in  the  capabil¬ 
ity  for  higher  dutch  roll  frequency  and  damping.  The  feedback  gains 
selected  for  each  configuration  are  also  shown. 


Effect  Of  Actuator  ttime  Constant  Changes 


VI.  Conclusions  and  Recommendations 


Conclusions 

It  has  been  seen  in  this  study  that  removal  of  the  vertical 
stabilizer  on  the  F-lll  results  in  an  unstable  aircraft.  The  analyses 
of  this  study  have  shown  that  lateral-directional  thrust  vector 
control  (TVC)  can  not  only  stabilize  the  F-lll  in  the  tail-off  con¬ 
figuration,  but  also,  can  achieve  a  better  dutch  roll  frequency  and 
damping  ratio  than  the  baseline  tail-on  aircraft .  For  the  two  flight 
conditions  examined,  the  response  of  the  tail-off  aircraft  to 
atmospheric  turbulence  was  as  good  as  or  better  than  the  tail-on 
aircraft.  The  price  that  must  be  paid  for  this  stability  is  a 
momentary  reduction  of  the  forward  thrust  due  to  the  deflection  of 
the  thrust  vector.  If  the  thrust  vector  is  allowed  higher  maximum 
excursions  and  rates,  the  dutch  roll  stability  is  better  than  that  of 
the  tail-on  aircraft.  The  problem  is  that  about  12$  of  the  forward 
thrust  (for  a  maximum  deflection  of  0.5  radians)  is  lost  momentarily 
at  the  full  deflection  of  the  thrust  vector.  It  was  seen  that  the 
maximum  TV  excursion  could  be  limited  to  as  low  as  n.2  radians 
(ll.5  degrees)  and  still  maintain  acceptable  flying  qualities.  This 
would  only  be  a  2%  loss  in  forward  thrust  for  full  deflection.  This 
assumes  a  worst  case  cruise  condition  (low  thrust)  in  an  environment 
represented  by  the  discrete  atmospheric  turbulence  model  at  the 
thunderstorm  levels.  But,  a  TV  excursion  rate  of  1.1  rad/sec 
(63  deg/sec)  is  required  with  the  above  excursion  value. 

Requirements  for  the  thrust  vector  deflection  and  rate  of 
deflection  were  specified  for  an  actuator  model  with  a  0.1  second 


time  constant.  It  was  seen  that  the  lowest  TV  rate  of  deflection  for 


acceptable  flying  qualities  was  1  rad/sec,  so  a  device  to  deflect  the 
thrust  vector  having  at  least  this  capability  would  be  required.  It 
was  seen  that  a  TVC  actuator  with  faster  dynamics  (than  the  0.1  sec. 
time  constant)  produced  a  better  dutch  roll  frequency  and  damping 
ratio.  The  slowest  actuator  possible  was  one  with  a  0.25  second 
time  constant.  Slower  actuator  dynamics  resulted  in  a  closed  loop 
coupled  roll-spiral  mode. 

The  control  system  feedback  gains  found  in  this  study  were  de¬ 
pendent  on  the  flight  condition  and  thrust  level,  implying  that  gain 
scheduling  is  necessary.  Since  it  was  shown  that  the  closed  loop 
system  was  very  sensitive  to  unmodelled  changes  in  feedback  gains , 
thrust  level,  and  actuator  dynamics,  the  gains  must  be  closely 
matched  with  the  actual  thrust  and  actuator  dynamics. 

Recommendations 

This  study  was  done  using  the  lateral-directional  small  pertur¬ 
bation  equations  of  motion.  It  is  recommended  that  the  study  be  con¬ 
tinued  using  the  non-linear  longitudinal  and  lateral-directional 
equations  of  motion.  Flight  conditions  such  as  a  landing  configuration 
and  maneuvering  flight  should  be  considered.  A  model  of  the  amount  of 
thrust  deflection  that  can  be  practically  realized  needs  to  also  be 
considered  in  further  analyses. 

A  problem  unique  to  utilizing  TVC  for  lateral-directional  stability 
is  one  of  an  engine  out  condition.  Since  the  F-lll  is  a  twin  engine 
aircraft,  this  would  help  minimize  the  problem,  but  there  still  may  be 
consideration  to  provide  a  small  vertical  surface  to  produce  a  more 
neutrally  stable  aircraft. 
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F-lll  TACT  Aircraft  Physical  Description 
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Appendix  2 

F-lll  TACT  Stability  Deriviatives 
Body  Axes 


MACH  0.8 

MACH  2.0 

CC 

=  90 

CC 

=  2° 

TAIL-OH 

TAIL-OFF 

TAIL-OH 

TAIL-OFF 

c 

*6 

(l/deg) 

-.0014 
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T» 
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cn 

(l/deg) 

.0012 

-.0012 

.0012 

-.0012 

c 

£p 

( 1/rad) 

-.30 

-.25 

-.11 

-.09 

c 

yp 

(l/rad) 

.07 

-.05 

.05 

-.04 

c 

n 

P 

(l/rad) 

-.03 

-.009 

-.01 

-.003 

C 

*r 

(l/rad) 

.31 

.29 

.017 

.086 

C 

yr 

(l/rad) 

0.0 

0.0 

0.0 

0.0 

cn 

r 

(l/rad) 

-.25 

-.08 

-.21 

-.07 

C 

\ 

(l/deg) 

.0001 

.0001 

.00013 

.00013 

C 

'‘a 

(l/deg) 

-.0008 

.  -.0008 

-.000067 

-.000067 

c 

n 

(l/deg) 

-.000067 

-.000067 

-.000067 

-.000067 

5 

a 
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MACH  0.8  MACH 

cc  =  9°  «  = 


TAILt-ON 

TAIL-OFF 

TAIL-ON 

*c 

\ 

( 1/deg) 

.003^ 

.001 

.0008 

\ 

(l/deg) 

.0003 

0.000 

.00005 

5v 

*c 

"a 

( 1/deg) 

-.0013 

-.0006 

-.00026 

*6  control  derivatives  in  Tail-on  Case 
r 


.0 

o 

TAIL-OFF 

.0007 

0.000 

-.000)4 
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Appendix  3 

Thrust  Vector  Control  Derivatives 

Thrust  vector  control  derivatives  are  required  in  the  same  form 
as  the  rudder  control  derivatives  in  order  to  directly  apply  them 
to  the  lateral  -  directional  small  perturbation  equations  of  motion. 

The  usual  sign  convention  for  rudder  control  is  that  left 
rudder  is  positive,  which  causes  a  negative  yawing  moment.  The 
same  sign  convention  will  be  used  for  the  thrust  vector  deflection 
in  that  a  thrust  deflection  that  causes  a  negative  yawing  moment 
is  considered  a  positive  angle  (6^).  This  is  shown  in  figure  3-1. 
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From  figure  3-1  "the  force  in  "the  Y  direction  due  to  thrust  "vectoring  is 

Y  =  T  sin  $ 
t  -v 

Using  a  small  angle  approximation: 

Sin  6=6  (rad) 
v  r 


Then: 


y  =i6 
t  v 


Since  the  engine  thrust  is  assumed  to  he  parallel  to  the  roll  axis 
in  the  X  -  Y  plane  the  rolling  moment  due  to  thrust  is  zero.  That  is: 

Lt=° 

The  yawing  moment  due  to  thrust  vectoring  is  (using  the  small  angle 


approximation) 


N  =  —  T  6  i 
t  -v  t 


■where  1^.  is  the  distance  between  the  aircraft's  center  of  gravity  and 
the  nozzle.  To  non-dimensionalize  the  thrust  vector  dimensional 
stability  derivatives, control  surface  derivatives  used  in  the  lateral- 
directional  small  perturbation  equations  are  defined  as  follows: 

y6  =  -a-s  V 


\  =  3, 1  *  C# 


Nr  q  s  b  C 

6=  uT  ns 


Now  let: 


Y '  =  T  6  =  m  Y.  6 
t  v  6  v 
v 


lt=° 


U.  =  -  T  6  1  =  -  Izz  N.  6 

t  v  t  6  v 

v 
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and  defining  the  thrust  vector  control  derivatives  as: 


\ 


v  m 
L 


6=0 

v 


N„  _  q  s  h  n* 

0  -  -  6 

V  Izz  V 

So  finally  the  non-dimensional  thrust  vector  control  derivatives 

C  m  Y6  _ 

yc  = - v  =  _T_ 

v  q  s  q  s 


=° 


C  IzzH6  T  1, 

n^ _ v  _  _ t 

v  q  s  1)  q  s  b 

For  steady  level  equilibrium  flight  conditions: 

T  =  D  =  CD  qs 


Then: 


Cy  Vs 

y6  =  — —  =  C„ 
,r  qs  D 


'*«  =0 


\  ,  Vs  H  _  _  CD  h 

v  qsb  b 


The  drag  coefficient  was  obtained  from  reference  8. 
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Listing  of  Integration,  ana  Plot  Program 
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Listing  of  Continuous  Random  Variance  Program 
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Unclassified. _  " 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (OTien  Data  Entered) 


DD  ,  55“,,  M73  EDITION  OF  I  NOV  65  IS  OBSOLETE 


UnrO  pqpj  f i 


SECURITY  CLASSIFICATION  OF  THIS  PACEfHTiwi  DtU  EnCr.d) 


he  eyalua,tion.  It  is  found  that  the  thrust  vector  control  produces 
a  response  as  good  as  or  better  than  the  baseline  aircraft. 

Requirements  for  the  thrust  vector  deflection  and  rate  of 
deflection  are  generated.  The  lowest  possible  rate  of  deflection  for 
acceptable  flying  qualities  is  shown  to  be  1  rad/sec.  er 


